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On  the  Computation  of  the  Eccentric  Anomaly,  Equation  of  the 
Centre  and  Radius  Vector  of  a  Planet,  in  Terms  of  the  Mean 
Anomaly  and  Eccentricity.  By  J.  Morrison,  M.D.,  M.A., 
Assistant  on  the  American  Ephemeris  and  Nautical  Almanac, 
Washington,  D.C. 

{Communicated  by  the  Foreign  Secretary.) 

Several  methods  of  computing  the  eccentric  anomaly  of  a 
planet  from  the  eccentricity  and  mean  anomaly  have  already 
been  given,  in  the  Monthly  Notices  and  elsewhere,  by  many  very 
able  mathematicians.  Some  of  these  methods  involve  consi-. 
derable  labour  even  after  the  formulae  have  been  deduced,  while 
others  give  only  rongh  approximations.  The  following  method 
of  treating  the  subject  is  due,  I  believe,  to  the  late  Professor 
Hansen,  as  indicated  by  him  in  the  Abhandlungen  der  Sdchs- 
ischen  Gesellschaft  der  Wissenfchafte>i,  Band  II.  It  gives  the 
eccentric  anomaly  with  great  accuracy  and  facility  by  means  of 
rapidly  converging  series,  and  in  most  cases  with  very  little 
labour,  especially  after  the  coefficients  of  the  several  terms  of  the 
series  have  been  computed  for  each  planet.  I  purpose  in  this 
short  paper  to  develop  the  subject  more  fully  than  is  done  in  the 
work  referred  to,  and  to  apply  the  results  to  the  computation  of 
the  eccentric  anomaly  in  the  case  of  each  of  the  primary  planets, 
as  well  as  of  the  equation  of  the  centre  and  the  radius  vector. 

Let 

M  =  the  Mean  Anomaly, 
E  =  the  Eccentric  Anomaly, 
V  =  the  True  Anomaly, 

e=the  Eccentricity, 
and 

€  » the  Napierian  base. 
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then  we  have  the  following  well-known  relations  : 

M  =  E-esinE;  (i) 

tan  Ai'=  .  A^ .  tan  lE.  (2) 

Hansen  assumes 

and 

that  is, 

y  =  cos  E  -f  V"^  sin  E, 
and 

2r  =  cos  M  +  V^T  sin  M, 

and  represents  the  coefficients  of  the  series  to  be  developed 
according  to  the  sine  and  cosine  of  multiples  of  E  by  P^^^  and 
Q'^\  thus 

y«=5    F^'^z"  (3) 

J=-0O 

^*  =  r"Q(?V  (4) 

'  A=-co 

where  h  and  i  are  integers,  and  it  is  evident  fi'om  the  form  of  the 
series  that  we  must  have 


nnA 

pT=p-T 

and 

Q^?=Q:f. 

except  i 

in  the  case 

of  h=o. 

Relation  between  P'^'  and  Q"i\ 

Put 

A  =  i  +  M  in  (3), 

and  we 

shall  have 

.V«'=2"    Ti%zi+>^. 

Multiply  by  z~''dM  and  integrate  between  the  limits  tt  and  —  ir, 
thus 


airv'-il  _,^=r 


y^-M-i-irf,  (5) 
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since 


Again,   in  (4)  pat  i=^  +  v,  multiply  hj  y~''dE,   and   integrate 
between  the  limits  w  and  —  tt,  thus 


^^^\     y-'-'^'^ 


..^^,^1,^---^^^         (^) 


since 

Let 
that  is 


i  +  ix=-k,  and  k  +  v  =  i; 


M  +  "  =  o. 
Then  (5)  and  (6)  become 


I       1«'^^ 

p^i'=— 7=^        y^-'-irf^,  (7) 


Integrating  (7)  by  parts,  we  have, 


_yi-'^z-'^dy,  (9) 


Comparing  (8)  and  (9),  we  have 

the  required  relation. 

Eliminating  E  and  M  from  the  equations 

y  =  c»V=r     2.  =  fMv'-T-  and  M  =  E-e  sin  E, 
we  get 
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which,  substituted  in  (8),  gives 

Put 

A  =  _  ^c/t, 

and  develop 

by  the  exponential  theorem,  thus 

.x(,-o=,+x(^-y-i)  +  ^^  (y-y-i)-^4-_^-(y-,-i)'+  etc. 

^        I-'        2^       6«       24=»       1202  ^   ••••••     ;y 

/.        \'  A.*  A'  A" 

1.2       2[3_      64^       4[6^  '* 

^  .  A''         A«          A«          \'    ,       A'»  .     , 

+  ( — • + +  — \  ifl 

\2_      [3_       2|4_      16_  ^24i6_  ^-^ 

<A'^         A<  A«  A»  A'»  . 

\L      ^       2[4,      I6_      24(1  '^ 

^  ,\^         A*     ,      A'  A«  A"  .     . 

+  ( ■ —   + +  --       .    • 1  V* 

^i3_      14.       2^      6\±      24[7_  ^^ 

[3.  "^14.       2[5.^6[6_      24|2_'^ ^^' 

+ 

+  ,^-(, ^+ ^! 

\m         I  .m+ 1       i.2.m+i.m  +  2 

+    .  .  .  )ir 


where  m  must  be  considered  a  positive  integer. 

Let  us  now  represent  the  coefficient  of  y"'  by  J^'"^  (Hansen's 
notation  for  the  Besselian  Function)  ;  then  the  above  series  may 
be  written  thus  : 

ni=oo 
«x(y-.v-'>  =  2        .T'TV".  (II) 

m=— » 

and  from  the  form  of  the  preceding  development  we  notice  that 
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From  the  last  two  equations  we  have, 

Q?  = 7=2        J(^M  y^-i-'-ldy.  (12) 

Now,  since 

when  /<.  is  an  integer  and 

I  y-lrfv  =  2r  v'  — I. 

we  must  have  h—i—m^o,  otherwise  (12)  would  vanish;  there- 
fore we  have 

2X-/-I    m=-oo       *     L-»-/^ 


hence 


But 


and  also 


=  2       jT; 


y  =  2        P'i'^=A  2         J'*I'V; 


y'  =  cosiE+  v^— isiniE 
5*  =  cos  AM  +  V^^  sin  AM  ; 

therefore  we  have 


A=» 


cos  iE+  >y— I  sin2E  =  i  2        J'  x '(cosAM+  v'— i  sin  AM). 

Equating  the  real  and  imaginary  parts  of  this  equation  we  have 

i  *=" 
cos  lE  -  -j  2        J^*!'")  cos  A  M  (13) 

siniE  =  ^2        J'*r"  sin  AM  (14) 


When  ^:=o,  the  second  members  of  the  last  two  equations  assume 
the  form  -;  hence  it  becomes  necessary  to  evaluate  them.  Puttinor 
/t=o  in  (7)  it  becomes 
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Differentiating 

we  get 

z-^dz=y~^dy —  {i+y~^)dy, 

which,  substituted  in  the  preceding  equation,  gives 

p«  =  -— :=  iyi-^  -  1 1-  I  yi-^)dy 

=  o, 

for  all  values  of  i  greater  than  i  or  less  than— i.  Therefore,  o,  i 
and  —  I  are  the  values  of  t,  which  make  the  second  members  of 

(13)  and  (14)  take  the  form  -. 
When  t=o,  we  have 

2ltV-l\_^^Z^  2         2 

y-^dy- — rr    -  (i  +y-*yy 


2itV-l\.^^-=X  2ir^-I     2   |j_„^^ 

=  1-0=1. 


When  1=1,  we  have 


When  i=:  —  I,  we  have 


^""■d^f  IT^,'*--- 1*--  ->-*)*=-»•■ 


For  all  other  values  of  i  we  have 

P^i^  =  o. 
From  (13)  we  find,  when  t=i,  2,  &c. 
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Reducing  the  coefficients  of  sin  M,  sin  2M,  &c.,  to  logarithms, 
we  have  the  following  series  for  each  of  the  primary  planets. 

Mercury, 
e  =  -20560478  fc"  =  42409''-03 


E  =  M  + [4-6251614]  s 
+  [2-8172458]  s 
+  [1-3733623]  s 
+  [0-0394355]  s 
+  [8-7587952]  si 
+  [7-5113252]  s 


p  M  +  [3-6333310]  sin  2  M 
n  3  M  +  [2-0747079]  sin  4  M 
n  5  M  +  [0-6972126]  sin  6  M 
n7M  +[9-3948531]  sinSM 
n9M  +[8-1318557]  sin  10  M 
n  1 1  M  + 


For  finding  the  true  anomaly  we  shall  have 

tan  ^  ;/=  [0-0905841]  tan  ^  E. 

Example.     Let  M=64°  10'. 

E  =  64°  10'  +  37969"-5526  +  337i"-9202- 142"-0962 

—  ii5"-5688—  i4''-9202  +  2"-i046  +  i"-095o 
+  o"-iii4— o"-o349— o"-oi32— o"-ooo8 
=  64°  io'+ii°  24'  32"-i5o  =  75°  34'  32"-i5o 
*y  =  87°  22'  2o''-29. 

Venus. 

e  = -00684331  .         e"=i4ii"-534 

E  =.M+  [3-1496888]  sin  M     +  [0-6839207]  sin  2  M 
+  [8-3942434]  sin  3  M  + [6-1783520]  sin  4  M 
tan  I )/  =  [0-0029720]  tan  |  E. 

The  Earth. 

«  = -01677 1 1  «"  =  3459"-287 

E  =  M  + [3-5389714]  sin  M     +[1-4624774]  sin  2  M 
+  [9-5620725]  sin  3  M  +  [7-7354525]  sin  4  M 
tan  ^  «/  =  [00072842]  tan  ^  E. 

Example.     Let  M= 7 1  °. 

•^  =  71° +  54'  3o"-7o6o+  i7"-8574-o"-i987-o"-0O527 
-71°  54'  48"-359. 
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Check. 

E=M  +  c"  sin  E 

=  7i°  +  3288"-359  =  7«°  54  48"-359- 

Mars. 
6=09326113  e"=i9236"-5 

E  =  M  +  [4-2836536]  sin  M  +  [2-9515369]  sin  2  M 
+  [17954344]  sin  3  M  +  [o-7i3o829]  sin  4  M 
+  [9-6715855]  sin  5  M  + [8  6560435]  sin6M 
+  [7-6585193]  sin  7  M  + 

tan  J  y  =  [00406208]  tan  i  E. 

Jupiter. 
e  = -0482519  e"  =  9952"-67 

E  =  M  +  [3-9978i3i]  sin  M     +[2-3800869]  sin  2  M 
+  [0-9384309]  sin  3  M  +  [9  5706569]  sin  4  M 
+  [8-2435235]  sin  5  M  +  [6-9424833]  sin  6  M 
tan  ^  »'=  [00209718]  tan  ^  E. 

Saturn. 

€  =  •0559428  e"=ii539"-02 

E  =  M  +  [4-0619992]  sin  M     +[2-5084302]  sin  2  M 
+  [1-1309242]  sin  3  M  +  [9-827i927]8in4M 
+  [8-5643061]  sin  5  M +  [7-3274187]  sin  6  M 
tan  ^  i/  =  [00243210]  tan  \ E. 

Urantts. 

•  =  •0463592  e"  =  9562"-2g 

E  =  M  +  [3-9804453]  sin  M     +  [2-3453568]  sin  2  M 

+  [08917040]  Sin  3  M+  [95011044]  sin  4  M 

+  [8-1567170]  sin  5  51  +  [6-8383189]  sin  6  M 

tan  ^  i/=  [0-0201480]  tan  ^  E. 

Neptune. 
e  = -0089903  e"=i854"'38 

E  =  M+  [3-2681948]  sin  M     +  [0-9209318]  sin  2  M 
+  [87497592]  sin  3  M  +  [66523734]  sin  4  M 
tan  ^  »'  =  [0-0039045]  tan  J  E. 
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Equation  of  the  Centre. 
From  the  theory  of  elliptic  motion  we  have 


dp       h 

dt       r'' 

where  h= 

Uv'.CJL- 

=hftlfthe 

area  described  in  a  unit  of  time ; 

T   _  2. 
irab       h' 

therefore 

whence 
therefore 

J      2rab , 

dv       Zirab 
dt  ~    Tr^ ' 

and 

(fM  ^  2ir 

dt  '^  T' 

Now, 

but 
and 
therefore 


dv    _  dv  ^    dt   _  ah 
dK'  ~di     dM.~  V-' 

r  =  a(i— e  cos  E); 
rfM       (i-ecosE)''' 


(15) 


Developing  the  second  member  by  the  binomial  theorem,  and 
changing  the  powers  of  cos  E  to  the  cosines  of  multiples  of  E, 
we  have 
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Greatest  Equation  of  the  Centre. 

When  the  equation  of  the  centre  has  its  maximum  value, 
d(v— M)=o  or  dv=dM;  therefore,  from  (15),  we  have 

cos  E  = i ^ 

e 

=  _L  +  3£!  +  7f!  +  77«'  + 
2-       2*        2'        2" 


Bin  E=l ;  - ^ ...  (10) 

2*        2"  2'* 

whence 

2         2*  2'^  2" 

„  0€*      4i7e*   .    io87e*  , 

sin  3  E  =  - 1  +21  +  tt_^  +  — ^  +    .  .  . 
J  2*         2"  2'* 

&c.,  &c. 

Any  one  of  these  will  give  the  eccentric  anomaly  when  the 
equation  of  the  centre  is  a  maximum.     Developing 

tan  Ji'=^i±f  .tanJE 
into  a  series,  we  have 


»'  =  E  +  2  tan  *  sin  E  + tan=-  sin  2E  +  -  tan  -  sin  3E+   &c., 
2  232 


where 


UBJ^  =  i-C-^)* 


=  -+  —  +  —  +  2—  +    &c. 
2         2'  2*  2' 

Substituting  the  values  of  sin  E,  sin  2E,  &c.,  tan*,  tan^X 
°  22 


&c.,  we  find 


.  =  E  +  .+  4^^    +iil^^   +^«'   +   ...  (17) 

2*.  3  2".  5  2'*.  7 
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Multiplying  the  second  series  of  (16)  by  e,  and  substituting  in  (i), 
we  get 

M  =  E-e+  —  +  49^5    +  H33g7   +   .  .  .  (ig) 

2.5       2"  2'* 

Subtracting  (18)  from  (17),  we  have 

=  2,+_ii_e3    +  J99.,s    +  17219,,   +   _  .  (,9). 

2*.  3  2'".  5  2'*.  7 

Reverting  this  series  we  have 


.C3 ^11— C' ^"^"-^     C^  -  .  .  .  (20) 

2       2«.3  2'«.3.5  22'.  3-^5.7 


Tlie  Radius  Vector. 

By  the  properties  of  the  ellipse  we  have  r  =  a  (i— e  cos  E), 
in  which  substitute  the  values  of  cos  E  ;  we  thus  deduce  the 
following   series,  a   check  upon  which  is  obtained  by  putting 

•  T 

M=o,  when  it  reduces  to  -~  =  i—e. 

a 


*  Delambre,  in  his  Astronomy,  vol.  ii.  p.  55,  has  given  this  series  correctly, 
but  farther  down  on  the  same  page  he  committed  an  error  in  reducing  the  term 

in  c',  which  he  gave  as    ',   "e'.     He  also  gave  the  second  term  of  series  (20) 
2    .  7 

incorrectly  as C.     Both  of  these  incorrect  values  were  copied  by  Baily 

2'.  3 
in  his  "  Astronomical  Tables  and  Formulae,"  and  from  him  by  several  other 
writers  on  Astronomy.  Mr.  Farrel  in  the  Monthly  Notices  of  the  Royal  Astro- 
nomical Society,  June  i860,  corrects  Baily  as  to  series  (20),  but  at  the  same 
time  adopts  Delambre's  incorrect  value  of  the  last  term  of  (19),  and  computes 
accordingly. 
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The  Logarithm  of  the  Radius  Vector. 

The   logarithm  of  the   radius  vector  may  be  expressed  in 
terms  of  the  eccentricity  and  mean  anomaly,  by  developing  the 

equation  -^i  —  e  cos  E  by  the  logarithmic  series,  thus. 


log  -=  -e  cos  E-—  cos*  E-—  cos»  E  -&c. 
a  2  3 


Transforming  this  series  from  powers  of  cos  E  to  the  cosines  of 
multiples  of  E  we  obtain  the  following  series,  which  expresses 

the  logarithm  of-  in  terms  of  the  eccentricity  and  the  eccentric 
anomaly. 
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Check,     Putting  M=o,  it  reduces  to 

r  e.''       c'       r*  f'" 

log-=- (/-+'+'+"+    .   .  .     +'^    +    .  .   .) 
a  234^  10  ' 

as  it  should. 

When  the  eccentricity  is  large,  as  in  the  case  of  tlie  orbits  of 
comets  and  some  of  the  asteroids,  the  preceding  series  for  the 
computation  of  the  eccentric  anomaly  are  not  sufficient!}'  con- 
verging'.    In  that  case  the  foUowing  method  may  be  used  : 

In  the  equation  M=E-e  sin  E,  let  E  =  M  +  .r,  that  is  x= 
e  sin  E  ;  then  we  shall  have 

M  =  M  +  .r-c  sin  (M  +  .r) 
or 

x  =  r  sin  ^l  pos  .r  +  e  sin  x  cos  M 

=  f  sin  M  (i 1-    .  .  .)  +  e  cos  M  (.r—  "-    f    .  .  .) 

2  6 

or  - 

■    TIT     /.  -ITS       f  sin  IVI  ,     e  00s  M  , 

e  sin  JVl  =  (i— e  cos  M).r4 —  .t'  +       -   -  .r^   —    .  .  . 

2  6 

Reverting  this  series,  we  have 

e  sin  E  =  .T 

_     e  sin  M        1/    ''  ■'''"  -^^    \'' 
I  —  c  cos  M     'Vi  —  «  cos  mJ 

therefore, 

E  =  M  +  e  sin  E 


—  M         P  sin  M        .  /    g  sin  M    V 
I  —  e  cos  M        \  I  —  e  cos  M/ 


It  will  generally  be  sufficient  to  use  only  the  first  two  terms  of 
this  series  which  will  give  the  first  approximate  value  of  E, 
which  we  will  denote  by  E'.  With  this  approximate  value 
compute 

M'  =  E'-.  sin  K'; 

then,  by  subtraction,  we  have 

M-M'  =  E-E'-c  (sin  E-sin  E') 

=  E  — E'  — e  (E  — ]*y)  cos  E',  approximately; 

whence 

M-M' 


E  =  E'  + 


l—e  cos  K' 


which  is  a  more  accurate  value  of  E. 

The  last  process  may  be  repeated  as  often  as  necessary,  but 
unless  the  first  approximation  is  very  far  from  the  truth  it  will 
not  be  necessary  to  proceed  with  the  computation  farther  than 
we  have  indicated. 
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This  is  the  same  as  the  method  recommended  by  the  late 
Professor  Encke  in  the  Berliner  Asironomisches  Jahrhuch,  1838, 
and  is  substantially  the  same  as  that  given  by  Grauss  in  his  Theoria 
Motus  Corporum  Coelestium,  art.  ii. 


Sff'ttinroodf  d-  Co.,  Prinltrs,  Nnr-Ureel  Square.  London. 
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THE  ORBIT  OF  THE  GREAT   COMET  ib)  1882. 

BY 

J.   MORRISON,   M.A.,   M.D., 

Assistant  on  the  American  Ephemeris  and  NatUical  Almanac,  Washington. 
{Communicated  by  the  Foreign  Secretary.) 

Erer  since  the  discovery  of  this  remarkable  comet  in  the 
early  part  of  September  18S2,  numerous  efforts  to  determine  the 
elements  of  its  orbit  have  been  made  by  astronomers  in  different 
parts  of  the  world,  and  in  no  similar  case  has  there  been,  we 
believe,  a  greater  disparity  in  the  results  obtained.  An  orbit 
like  the  one  under  consideration,  which  is  so  situated  with  re- 
spect to  the  orbits  of  the  planets  that  the  effect  of  planetary 
perturbations  is  absolutely  inappreciable,  must  furnish  a  good 
example  of  pure  and  undisturbed  elliptic  motion,  and  therefore 
we  would  naturally  expect  that  a  system  of  elements  computed 
from  three  or  four  good  observations  suitably  selected  would 
very  nearly,  if  not  exactly,  represent  its  motion.  In  this  ex- 
pectation, however,  we  are  disappointed,  for  no  orbit  has  hitherto 
been  published  which  completely  satisfies  the  observations,  espe- 
cially those  made  after  the  middle  of  October  or  thereabouts. 
From   the  time  of  discovery  until  about  this  date  the  nucleus 
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presented  a  tolerably  well-defined  circular  disk  which  could  be 
as  easily  observed  as  that  of  a  planetoid  ;  but  soon  afterwards 
extraordinary  changes  began  to  take  place  in  the  structure  of  the 
nucleus,  which,  from  being  circular  in  outline,  became  oval  or 
elliptical,  with  from  two  to  four  or  five  bright  points  or  centres 
of  condensation.  Thus,  from  about  October  20  until  May  26 
the  length  of  the  nucleus  varied  from  57"  to  about  135",  and 
the  position-angle  of  its  major  axis  varied  from  113°  to  255° 
— which  is  nearly  at  right  angles  to  the  meridian — hence 
it  is  not  difficult  to  see  that  an  error  of  30",  or  even  more,  in 
R.A..  may  exist  in  any  observation  made  about  this  time,  accord- 
ing to  the  point  of  the  nucleus  observed.  The  error  in  Declina- 
tion will  he  much  less  than  that  in  R.A.  for  the  reason  given 
above.  The  following  two  observations  made  at  the  Litchfield 
Observatory,  Clinton,  N.Y.,  and  at  the  Naval  Observatory,  Wash- 
ington, are  here  adduced  to  show  the  inconsistency  of  the  results 
obtained  when  difierent  points  of  the  nucleus  were  observed  : — 


■Washington  Mean  Time. 

R.A. 

Dec. 

1882,  Oct.    24712191 

h     m        s 
10     6  2578 

- 17    6  482 

24739461 

10     5  52*29 

-17     5     6-8 

According  to  the  geocentric  motion  of  the  comet  at  that  time 
these  observations  are  inconsistent  with  each  other,  and  there- 
fore of  little  or  no  value  in  the  present  investigation.  That  the 
discrepancy  is  due  to  this  cause  alone  is  evident  from  the  fact 
there  were  no  errors  in  the  position  of  the  comparison  star  or  in 
the  reduction  of  the  observations,  which  were  re-examined  by  the 
observers.  It  is  therefore  not  at  all  suprising  that  periods 
differing  very  widely  from  each  other  have  been  obtained  by 
different  computers,  according  to  the  observations  employed. 
Thus,  from  post-perihelion  observations  Mr.  Chandler  obtained 
periods  of  31 15  days  and  4070  years  respectively,  and  between 
these  extremes  we  have  periods  of  269,  652,  712,  793,  843,  997, 
and  1376  years.  The  most  probable  period  is  between  712  and 
793  years,  as  will  appear  presently.  With  the  view  of  com- 
paring the  elements  obtained  from  observations  made  aft^er  peri- 
helion passage  with  those  deduced  from  observations  made  before 
and  after  that  epoch,  we  select  the  following  observations  : 


Washington 

Mean  Time. 

B.A. 

Dec. 

Where  made. 

I 

1882,  Sept. 

d 

8  1 13405 

144  16  52"-5 

0      /        // 
-  0  57  464 

Windsor 

2 

8452672 

144  59  514 

-  0  56  30  0 

The  Cape 

3* 

14-348690 

161  28  200 

+  04  129 

Melbourne 

4* 

i6'368673 

170  39  263 

+    I  12    48 

,, 

5 

16-936700 

Ingress  on 

Sun's  disk. 

The  Cape 

6* 

19969786 

168  34  57-6 

-  0  34  285 

W.ijihingfon 

7 

Oct. 

8720457 

>57     I  37-9 

—  lo  40  219 

ji 
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Washington  Mean  Time. 

K-A. 

Dec. 

Where  made. 

d 
8                        27710103 

150  31    33-9 

-°i8  14  311 

Clinton 

g*              Nor.   1 5740418 

141    57  408 

-24  49  189 

Washington 

10*                        24700923 

136  34    33 

—  27    21    267 

„ 

II*              Dec.    11619921 

124    5  253 

—  30    16    281 

>• 

12       1883.  Jan.    13778694 

100    3  549 

-26  48   395 

Windsor 

13                 Feb.  28304652 

87  56  iSo 

-15  16    7-1 

Washington 

14                 Apr.     6097128 

89  31  28  9 

-94  49-2 

Palermo 

15                 May     6-252250 

93  32  126 

-  6  21  303 

Cordoba 

16                           26257880 

97  II  126 

-  5  30    00 

.. 

These   observations 

embrace    nea 

,rly 

the    entii-e 

period   of 

visibility  of  the  comet. 

No.    7    was 

made  with  the 

equatorial, 

the  comparison  star  being  a  well-known  one,  whose  position  was 
subsequently  verified  by  the  Transit  Circle;  Xo.  8  is  marked 
by  the  observer  as  "  very  good  ;  "  in  Nos.  g  and  1 1  the  "  main 
point  of  condensation  near  the  following  end  of  the  nucleus  was 
observed,"  and  in  13  "  the  middle  point  between  the  second 
and  third  points  of  the  nucleus."  Xo  remarks  respecting  the 
remaining  observations  are  made  by  the  resp>ective  observers. 
From  Nos.  6,  7,  and  13  we  have  obtained  the  following 
system  of  elements  : — 

I. 
T         1882,  Sept.  1700858  Washington  M.  T. 

M  69  37  277  Mean  Equinox  of  1882-0, 

Q  346     I  38  8 

«  38    o     7-7 

log.  q  78922566 

log.  e  9-9999660 

log.  a  19992346 

P  997'36  years.     Motion  retrograde. 

These  elements,  although  computed  with  great  care,  do  not 
represent  the  motion  in  the  neighbourhood  of  perihelion.  The 
residuals  for  the  middle  place  are  much  greater  than  what  would 
likely  be  due  to  the  unavoidable  errors  of  computation  and  obser- 
vation. Among  quite  a  number  of  places  computed  the  smallest 
residuals  obtained  were  the  following  : — ■ 


AACOS/S 

Afl 

3 

Sept  14-34  + 

u 
+  27-21 

+  17-20 

9 

Nov.  15-74  + 

-2449 

.  +12-53 

It  is  very  probable  that  a  slight  error  in  the  third  place  has 
to  some  extent  vitiated  these  results.  The  interval  between  the 
first  and  third  observations  is  161-3247  days,  during  which  time 
the  comet  described  an  arc  of  only  16°  i2'-i,  an  amount  of 
motion  far  too  small  to  permit  of  an  accurate  determination  of 
*  Meridian  observations. 
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the  tleraents,  especially  in  orbits  like  the  present  one,  whose 
eccentricity  is  very  large. 

For  the  computation  of  a  second  system  of  elements,  we  will 
take  observations  which  give  a  large  amount  of  motion,  and 
therefore  select  the  three  meridian  observations  of  September  14, 
September  19,  and  November  15;  the  last  being  the  first 
meridian  observation  made  after  the  comet  had  ceased  to  be 
visible  in  daylight,  and  had  passed  from  the  immediate  vicinity 
of  the  Sun.  At  this  time,  too,  the  process  of  disintegration  in 
the  nucleus  had  not  advanced  very  far,  and  from  the  remark 
appended  to  the  observation  it  is  very  probable  that  it  is  very 
nearly  as  accurate  as  either  of  the  other  two  about  which  there 
is  no  doubt. 

In  order  to  show  the  degree  of  accuracy  to  which  the  com- 
putation has  been  carried  we  here  give  the  results  of  the  last 
two  hypotheses.  The  notation  is  the  same  as  that  used  by 
Gauss  in  his  Theoria  3Iotus,  with  the  exception  that  we  have 
employed  u  to  denote  the  argument  of  latitude  and  p  the  curtate 
distance.  The  ratio  of  the  sector  to  the  triangle  corresponding 
to  the  first  and  second  observations,  and  also  to  the  first  and 
third,  is  negative  and  numerically  less  than  unity. 


1. 

IL 

log.P 

9'5492503'» 

9-5492433« 

log.  Q 

8-44584o8» 

8-4458447» 

01  +  «r 

2  25  49788 

2  25  5o'oi6 

log.  Qc  sin  w 

0-3691234 

0-3691230 

z 

48  i  i'^43 

0  /   // 
48  I  I  -936 

log.  r' 

9-3402558 

9-3402558 

log.  p' 

00539758 

00539758 

log.  p 

99972618 

99972631 

log.  p" 

0-0734265 

0-0734271 

log.  r 

93097846 

9-3097852 

log.r" 

0-2226261 

0-2226265 

Si 

346  0  4r38s 

346  0  4i"-254 

i 

38  0  7-03 

38  0  7-19 

Uu"-u) 

164  52  45-703 

164  52  47-13 

iK-«) 

157  55  32977 

157  55  3427 

i(M"-tt') 

6  57  12-726 

6  57  12-86 

log.  V 

0-1412227 

0-1412226 

log.  1," 

9-5876330M 

958763i9» 

log.  p. 

9-5492420» 

9'5492433» 

log.  q; 

8-4458464n 

8-445845on 

a: 

+  0-0000083 

00000000 

y 

-  00000056 

-  0  0000003 
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The  three  values  of  log  p,  the  semi-parameter,  are  each  equal 
to  8"  1 905 1 56,  which  proves  the  accuracy  of  the  entire  computa- 
tion up  to  this  point.  On  completing  the  computation  of  the 
elements  from  these  results,  we  have  the  following  system : — 

n. 

T  1882,  Sept.  17010163  Washington  M.T. 

M  69  35  30483  Mean  Equinox  of  1882*0. 

Q  346     O  41  254 

«  38     o     719 

log.  q  78895067 

log.  0  9  9999578 

log.  a  1-9016841 

P  7i2"io  years.     Motion  retrograde. 

The  residuals  for  the  middle  place  are  zero.  It  now  remains  to 
be  seen  how  closely  these  elements  satisfy  the  above  observa- 
tions. 

The  residuals  are  given  in  the  following  table,  where  C — O 
denotes,  as  usual,  the  difference  between  the  compxUe'l  and  the 
observed  places,  and  A  the  distance  from  the  earth. 


C-0. 

dAcos/S 

dp 

log.  A 

I 

-    371 

-   0-28 

00604274 

2 

-    075 

-p     I -02 

0-0569198 

4 

+    i-6o 

+    082 

9-9900265 

7 

—  016 

-    044 

01406414 

8 

+    i-oi 

+    037 

0-1662759 

10 

+  092 

—    1-12 

0-1766315 

II 

+  13-41 

+   6'45 

01916736 

12 

—  6120 

+  3752 

0-2798283 

•3 

-45-68 

+   382 

0-4629200 

14 

—  5212 

-2724 

0-5880931 

»5 

-6407 

-31-20 

0-6640691 

16 

-5412 

-2074 

0-7022010 

Traiisit  over  the  Suit's  Disk. 

The  ingress  on  the  Sun's  disk  took  place  on  Septembei 
17^  4**  50™  58',  Cape  mean  time,  according  to  Mr.  Finlay,  and  at 
^h  ^Qin  22»,  according  to  Dr.  Elkin  {Montidy  Notices,  November 
1882).  Adopting  the  former  and  applying  the  correction  for 
aberration,  which  is  —  o'*-oo5698,  we  have  for  the  position  of  the 
comet  the  following  results  : — 
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Washington  Mean  Time.  A  0  log.  A 

1882,  Sept.  16-931001         174  27  1983         -o  13   12  39  99952130 
Sun's  position,  ©    ...         174  36  1906         — o    o    0'22 


Difference 


o     8  59-23 


13  12-17 


Distance  between  the  centres 
Sun's  semi-diaraeter 

Comet's  distance  from  Sun's  limb 


15  58-28 
IS  5773 

055 


The  diagram  shows  the  position  of  the  comet  with  respect  to 
the  San  at  the  time  of  ingress. 

N  S  is  a  meridian ;  C  D  the  apparent  path  of  the  comet ; 
N  AC  the  angle  of  position,  122°  2^'"] 
from  the  north  point  of  the  Sun's  limb 
towards  the  west. 

Since  this  second  system  of  elements 
represents  the  motion  of  the  comet  with 
great  accuracy  from  the  time  of  dis- 
covery till  about  the  middle  of  Novem- 
ber, and  from  this  date  until  the  last 
observation  here  recorded,  with  all  the 
precision  that  can  be  expected  under 
the  circumstances,  we  may  also  give  the 
following  facts  as  tending  to  show  more 
strongly  perhaps  the  characteristics  of  this  remarkable  body. 

At  the  time  of  ingress  on  the  Sun's  disk  the  comet's  true 
anomaly  was  92"  7'  29"'668,  which  was  described  in  one  hour  and 
fifty-four  mimites.  From  the  time  of  the  first  observation  made 
at  Windsor  by  Mr.  Tebbutt  until  May  26th,  an  interval  of 
260-1447  days,  the  comet  described  340°  22''7  of  its  orbit,  but  it 
will  require  more  than  seven  centuries  to  describe  the  remaining 
19°  37''3-  Taking  the  mean  solar  parallax  at  8'''848,  the  value 
adopted  in  the  English  and  American  Nautical  Almanacs,  the 
comet's  perihelion  distance  was  716200  miles  from  the  Sun's 
centre,  or  285500  miles  from  the  surface. 

The  velocity  at  perihelion  was  295-56  miles  per  second,  which 
is  only  -0160  miles  per  second  less  than  that  due  to  a  parabola 
at  the  same  perihelion  distance. 

When  the  comet  reaches  aphelion,  it  will  have  plunged  down 
into  space  a  distance  equal  to  1 60  radii  of  the  Earth's  orbit,  or  5^ 
times  the  radius  of  Neptune^s  orbit,  at  which  enormous  distance 
its  velocity  will  be  only  75  feet  per  second. 

At  the  time  of  Dr.  Elkin's  first  observation  on  Sept.  8,  the 
apparent  diameter  of  the  coma  was  estimated  at  from  40"  to  50", 
and  the  strongly  condensed  cone  or  nucleus  in  the  centre  at 
from  10"  to  15".     Assuming  the  means  of  these  measurements, 
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we  easily  find  the  diameter  of  the  coma  at  that  time  to  have 
been  22500  miles,  and  of  the  nucleus  7600  miles,  or  a  little 
less  than  that  of  the  Earth.  As  the  comet  approached  the  San, 
the  nucleus  contracted,  and  about  half  an  hour  before  its  dis- 
appearance on  the  solar  disk  Mr.  Finlay  found  by  two  micro- 
metric  measures  the  apparent  diameter  to  be  4",  while  Dr.  Elkin 
estimated  it  at  some  5'.  The  former  of  these  gives  a  diameter 
of  1770  miles,  and  the  latter  2215  miles,  or  a  little  greater  than 
that  of  the  Moon. 

The  second  system  of  elements  here  presented  proves  con- 
clusively that  an  orbit  can  be  found  which  satisfies  the  motion 
of  the  body  equally  well  both  before  and  after  perihelion  ;  and 
therefore  we  may  safely  conclude  that  no  ret;\rdation  occurred 
during  the  perihelion  passage;  or,  if  there  wis,  it  is  absolutely 
inappreciable.  There  are  two  other  systems  of  elements  which 
diifer  considerably  in  the  periodic  time  from  oar  second  system, 
but  the  agreement  of  some  of  the  other  elements  is  very  re- 
markable. We  refer  to  those  of  Prof.  Frisby  and  Dr.  Kreutz ; 
the  former  founded  his  calculations  on  the  observations  of 
Sept.  19,  Oct.  8,  and  Xov.  24,  and  the  latter  on  normal  places 
derived  from  observations  extending  from  S-^pt.  8  to  Nov.  15. 
The  three  systems  just  referred  to  are  hert'  given  for  com- 
parison. 
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log.  e 

log.  a 
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Prof.  Frisby. 
1882,  .Sept.  170088 

6q  36  1279 

346     I     791 

38     o     7-84 

78904739 

99999606 

I  9331366 

793869  yrs. 

.  Mean  Equinox  18820   Mean  Equinox  1882-0    Mean  Equinox  1882*0 
*  Washington  Mean  Time.         +  Motion  retrograde 


Dr.  Morrison. 
18S2,  Sept.  17010163 

^  35  3o''48 

346    o  41-25 

38    o    719 

7-8895067 
99999578 
1-9016841 

712-10  JT8. 


Dr.  Kreuta. 

1882,  Sept.  17-00993 

o      /       ^ 
69  36     150 

346     I   27-20 

38    o  19-9 

7-8894760 
9-9999610 
19505596 
843-10  yrs. 


This  comet  is  probably  identical  with  the  comet  of  370  B.C., 
as  already  pointed  out  by  Mr.  Maxwell  Hall.  If  we  divide  the 
interval  between  370  B.C.  and  1882  a.d.  by  three  we  get  a  period 
of  751  years,  which  is  very  nearly  the  arithmetic  mean  between 
our  period  and  that  of  Prof.  Frisby.  This  period  would  give  a 
return  in  a.d.  1131  or  1132,  in  both  of  which  years  a  great  comet 
appeared  (^Chinese  Annals). 
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The  Orbit  of  Pons'  Comet.  By  J.  Morrison,  M.A.,  M.D.,  Assistant 
on  the  American  EpJiemeris  and  Nautical  Almanac,  Washing- 
ton, D.C. 

The  only  orbits  of  this  comet  which  have  been  pnblished  are 
those  of  the  late  Prof.  Encke  and  of  Messrs.  Schulhof  and 
Bossert :  the  former  deduced  from  observations  extending  over 
a  period  of  only  two  months  in  the  year  1812,  and  the  latter 
from  observations  made  in  1883,  and  subsequently  corrected  two 
or  three  times  by  equations  of  condition.  Neither  Encke's 
orbit  nor  that  of  Schiilhof  and  Bossert  represents  the  motion  of 
the  comet  with  that  degree  of  accuracy  which  modern  investiga- 
tions demand.  In  order  to  obtain  an  independent  determina- 
tion of  the  elements,  I  select  the  three  following  observations 
made  at  Washington : — 

Washington  M.T.  Apparent  a  Apparent  5 

1883  Oct.  10-307425      248  25  56-1      +5^  51  49-0 
Dec.  27-249330      317  34  504      +29  31  14-8 

1884  Jan.  21-299233      357  36  IS  -13  31  23-1 

The  comparison  stars  were  respectively,  a  Berliner  JaJirbuch 
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star ;  t,  Gygni ;  and  a  yth-magnitTide  star  in  Cetiis,  as  given  in 
the  Catalogues  of  Weisse,  Lamont,  and  Stone. 

From  these  observations  I  obtain  the  following  system  of 
elements ; — 

T     1884  Jan.  25729107,  G.M.T. 

•     199  12  5048   I  Mean  Equinox 


Q 

254    9  4582   )    0 

( 

74    2    073 

e 

-95414506 

logs 

9-8896779 

log  a 

1-2282918 

P 

69-572  years.     Motion  Direct. 

The  residuals  for  the  middle  place  are  zero.  It  will  be 
noticed  that  these  elements  diflFer  considerably  from  those  of 
Encke  and  Messrs.  Schulhof  and  Bossert.  The  following 
residuals  show  how  closely  this  system  of  elements  satisfies  the 
observations  made  from  Sept.  5,  1883,  to  Feb.  5,  1884: 


Date. 
Washington  iI.T. 

C- 

COBfidX 

-0 

d» 

log  A 

Obeerration 
made  at 

[883  Sept. 

511  + 

-  3-57 

-372 

0-3668189 

Milan. 

23-32  -H 

-   365 

-5-09 

O-32S2SS3 

Washington, 

Oct. 

31-29  + 

—  6-14 

-6-59 

0-2103271 

„ 

Nov. 

133  + 

+  373 

—  6-90 

0-2060407 

i> 

27-29  + 

+  4-82 

+  7-90 

0  0707879 

„ 

Dec. 

6-25  + 

+ 10-03 

-I  18 

00093336 

„ 

2124  + 

+ 10-60 

-720 

98932833 

„ 

[884  Jan. 

3-23  + 

-  312 

—  2-29 

9-8126678 

„ 

9-30  + 

—  0-26 

-1-48 

98023359 

,, 

17-30  + 

-  0-82 

-4'30 

9-8246773 

„ 

25-29  + 

-    I-I2 

+  i-io 

9-8749017 

„ 

Feb. 

1-28  + 

+    1-78 

+  0-14 

9-9258903 

»» 

5-01  + 

-    274 

+  1-23 

9-9522952 

Rome. 

Owing  to  errors  in  the  adopted  position  of  the  comparison 
stars,  many  of  which  have  not  been  very  accurately  determined, 
and  also  to  the  uncertainty  always  existing  in  observations  of 
comets  in  consequence  of  the  ever  varying  form  of  the  nucleus, 
these  residuals  are  less  than  what  we  may  reasonably  expect  to 
arise  from  these  causes  alone,  not  to  mention  the  unavoidable 
errors  of  computation,  even  when  seven -figure  tables  are  em- 
ployed, as  in  the  present  instance.     The  formation  of  equations 
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of  condition  for  the  purpose  of  correcting  these  elements,  cannot 
be  here  applied  with  any  advantage. 

The  only  planets  which  can  have  any  appreciable  influence 
on  the  motion  of  the  comet  during  the  time  it  was  in  the 
neighbourhood  of  perihelion  are  Venus  and  the  Earth.  Assuming 
the  elements  I  have  obtained  above  as  the  osculating  elements 
for  the  epoch  Jan.  26*0,  Greenwich  Mean  Time,  the  perturba- 
tions of  the  rectangular  co-ordinates,  arising  from  the  action  of 
Venus,  are  found  to  be  as  follows  : — 


Greenwicli  M.T. 

&x 

«y 

s« 

1884  Jan. 

26-0 

—  o-o 

-  o-o 

+  0'0 

31-0 

— o-i 

-  0-3 

+  01 

Feb, 

5-0 

—  0"2 

—   i-o 

+  0-4 

lO'O 

-0-4 

-    2-0 

+  09 

15-0 

-06 

-  3-2 

+  17 

20*0 

-09 

-  47 

+  2-6 

25-0 

-13 

-  6  5 

+  3-8 

Mar. 

i-o 

-17 

-  8-5 

+  50 

60 

—  2-1 

-IO-8 

+  6-5 

II-O 

-2-6 

-13s 

+  80 

which  are  expressed  in  units  of  the  seventh  decimal  place.  These 
perturbations  do  not  sensibly  affect  the  inclination,  and  their 
influence  on  the  node  is  less  than  a  second  of  arc  ;  they  are  then 
almost  inappreciable,  and,  indeed,  they  scarcely  exceed  the 
uncertainty  of  the  seventh  decimal  place  in  the  final  results  of 
extended  computations  such  as  the  one  now  under  considera- 
tion. The  perturbations  of  the  rectangular  co-ordinates  arising 
from  the  action  of  the  Earth  are  even  less  than  those  produced 
by  Venus,  by  reiason  of  the  Earth  being  at  a  greater  distance 
from  the  comet. 

According  to  my  elements,  the  aphelion  distance  is  33'056, 
which  exceeds  Neptune's  mean  distance  from  the  Sun  by  about 
three  mean  radii  of  the  Earth's  orbit,  and  although  the  inclina- 
tion is  large  (74°  2'),  the  comet  will  nevertheless  be  sensibly 
affected  by  the  action  of  all  the  superior  planets  during  the 
greater  part  of  its  period.  The  amount  of  this  action  during 
the  last  revolution  may,  perhaps,  be  approximately  found  by 
comparing  my  elements  with  those  of  Encke,  as  given  in  the 
Annales  de  I'Ohservatoire  de  Bruxelles,  tom.  i.  If  we  refer 
Encke's  elements  to  the  mean  equinox  and  ecliptic  of  i884'o, 
and  place  them  beside  mine  in  order  to  facilitate  comparison,  we 
have  the  following  results  : — 


884. 

of  Pons' 

Comet. 

Morrison. 

Encke. 

T* 

1884  Jan.  2572917 

1812,  Sept  15-313553 
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199  12  50-48 

0      1       II 
199  17     589 
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254     I  28-91 

i 

74     2     073 
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logq 
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69572  years 
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Assuming  these  valaes  of  T  to  be  correct,  the  actnal  period  of  the 
last  revolution  was  71  "3605  years,  and  therefore  Encke's  period 
has  been  increased  by  247-09  days.  It  is  very  probable,  how- 
ever, that  Encke's  elements  are  considerably  in  error,  for 
Messrs.  Schnlhof  and  Bossert  have  deduced  from  the  original 
observations  of  181 2  the  following  elements,  which  agree  tolerably 
well  with  mine  : — 


»—  ii=eo 


e 
logq 


I      II 
12  32-5]  Jlean  Equinox, 
6  15-^  j  18800 


199 

254     6  15-3 
74    3  20-4 
0-9549960 
98893650 


If  we  add  3'  2o''*92,  or  four  years'  precession,  to  S8,  we  get 
254°  9'  36"'2,  which  is  almost  identical  vrith  my  value,  while  the 
values  of  w  and  i  will  not  be  sensibly  affected  by  this  change. 


♦  Ctreenwich  Mean  Time. 
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The  Apparent  Orbit  of  a  Satellite  of  a  Superior  Planet.  By 
J.  Morrison,  M.D.,  M.A.,  Ph.D.,  Assistant  on  the  American 
Ephemeris,  Washington,  DC. 

The  apparent  orbit  of  the  satellite  of  a  superior  planet  is  the 
pi'ojection  of  the  real  orbit  on  a  plane  perpendicular  to  the  line 
joining  the  Earth  and  planet,  and  is,  therefore,  in  general  an 
ellipse.  In  this  paper  I  shall  assume  that  the  elements  of  the 
orbit  of  the  satellite  are  known,  and  proceed  to  develop  as 
much  of  the  subject  as  is  necessary  for  the  computation  of  an 
ephemeris. 

Let 

a,  5,  A     be  the  polar  coordinates  of  the  planet  as  seen  from  the  Earth  ; 
a',  5',  A'  the  polar  coordinates  of  the  satellite  n.s  seen  from  the  Earth ; 
«,  d,  r     the  polar  coordinates  of  tlio  satellite  as  seen  from  the  planet, 

referred  to  the  equfitor  and  equinox  of  the  epoch  for  which  the 
ephemeris  is  to  be  computed. 

The  general  expressions  for  the  rectangular  coordinates  in 
tSfms  of  polar  coordinates  are 

X  =  A  cos  5  cos  o 
If  =  A  COS  5  sin  a 
£  —  A  siu  S. 
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Therefore  by  equating  the  rectangular  coordinates  we  hare 
a'  cos  8'  cos  a'  =  A  cos  5  cos  a+r  cos  d  cos  a 
A' cos  8*  sin  a*  =  A  cos  5  sin  a  +  r  cos  (f  sin  a  ...         (i) 

A'  sin  8*  =  A  sin  8  +  r  sin  d. 

Multiplying  the  first  of  this  group  by  cos  a,  the  second  by  sin  a, 
and  adding  the  results,  then  multiplying  the  first  by  sin  a,  the 
second  by  cos  a,  and  subti-acting,  we  have 

A'  cos  8*  cos  (a'  — a)  =  A  cos  S  +  r  cos  d  cos  (a— a) 

A' cos  8*  sin  (a' —  a)  =  r  cos  d  sin  (a  —  o)  ...         ...         (2) 

A'  sin  8*  =  A  sin  8  +  r  sin  d. 

We  have  now  to  consider  the  spherical  triangle  formed  by 
the  north  pole  of  the  equator,  the  planet,  and  the  satellite  ;  thus 
let  Q  be  the  north  pole,  P  the  planet,  and  S  the  satellite.  In  the 
spherical  triangle  Q  P  S  we  have 


PQ  =  90"  -  8,  QS  =  90°  -  8'; 
and 

PQS=a'-a. 

Put 

PS  =  s,  the  apparent  angular  distance  of  the  satellite  from  the  planet 
as  seen  from  the  E<irth,  and 

QPS  —  p,  the  angle  of  position  of  the  satellite  reckoned  from  the  north 
toward  the  east  and  from  0°  to  360°. 

By  the  fundamental  formulae  of  spherical  trigonometry  we  have 
sin  s  sin  p  =  cos  8'  sin  (o'  —  o) 

sin  s  cos  p  =  cos  8  sin  8'  — sin  8  cos  8*  cos  (^a'—a)  ...         (3) 

cos  8  =  sin  8  sin  8*  +  cos  8  cos  8'  cos  (a'  —  a) 

Multiplying  by  A'  and  substituting  from  (2)  this  group  becomes 

A'  sin  s  sin  p  =  r  cos  d  sin  {a  — a) 

A'  sin  s  cos  p  =  r  {cos  8  sin  d  — sin  8  cos  d  cos  (a  — a)]     ...         (4) 
A'  cos  «  =  A  +  r  {sin  8  sin  d  +  cos  8  cos  d  cos  (a— a)], 

Arom  which  we  must  eliminate  a,  d,  and  A'. 
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Lot  V  D  represent  an  arc  of  the  equator,  N  S  an  arc  of  the 
orbit  of  the  satellite,  E  the  Earth,  P  the  phmet,  and  S  the 
satellite  ;  draw  the  arc  S  D  perpendicular  to  V  D,  and  let  E  V 
be  the  direction  of  the  vernal   ec^uinox.     Then   N  will   be  the 


position  of  the  ascending  node  of  the  orbit  of  the  satellite  on  the 
equator  ;  V  D  and  S  D  the  planetocentric  Right  Ascension  and 
Declination  of  the  satellite,  which  we  have  denoted  by  a  and  d 
resppctivelj. 
Put 

VN  =  N,  the  right  ascension  of  the  ascending  node ; 

SND  =  i,  the  inclination  of  the  orbit  to  the  plane  of  the  equator ; 

NS  =  ti,  the  angular  distance  of  the  satellite  from  the  ascending  node ; 
then 

ND  =  a-N. 

From  the  right-angled  spherical  triangle  N  S  D,  whose  centre  is 
P,  we  have 

cos  d  cos  (a  — N)  —  cos  u 

cos  £?  sin  (a  — N)  =  sin  «  cos  ?"  ...         ...         ...         (5) 

sin  d  =  sin  n  sin  i. 

Multiplying  the  first  of  this  group  by  cos  (a  — N),  the  second  by 
sin  (a  — N),  and  adding  the  results  ;  also  multiplying  the  first  by 
sin  (a— N)  and  the  second  by  cos  (a— N),  and  subtracting,  we 
have 

cos  d  cos  {a  — a)  =  sin  u  cos  »  sin  (a— N)  +  cos  u  cos  (o— N) 

cos  d  sin  {a  — a)  =  sin  u  cos  i  cos  (o  — N)  — cos  u  sin  (o  — N) 
sin  d  —  sin  «  sin  i, 

which  substituted  in  (4)  give,  after  simple  reductions — 
A'  sin  s  sin  p  =  r  {sin  u  cos  i  cos  (a  — N)  —  cos  11  sin  (a  — N)} 

A'  sin  s  cos  p  =  r  {sin  u  [cos  5  sin  i  — sin  S  cos  i  sin  (a  — N)] 

—  sin  S  cos  u  cos  (o— N)} 

A'  cos  s  =  A  +  r  {sin  u  [sin  5  .sin  t'  +  cos  8  cos  /  sin  (o— N)] 

+  cos  5  cos  u  cos  (a  —  N) } , 
from  which  A'  must  be  eliminated. 
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In  the  plane  triangle  EPS  formed  by  joining  the  Earth, 
planet,  and  satellite,  let  a-  denote  the  sapplement  of  the  angle 
EPS;  then  from  the  diagram  (fig.  2 )  we  have 

A'  sin  s  =  r  sin  ff 
and 

a'  cos  6  =  r  cos  (T  +  A        ...         ...         ...         (6) 

which  being  substituted  in  the  equations  of  the  last  group,  give 
after  reduction 

sin  <r  sin  p  =  sin  u  cos  1  cos  (a  — N)  — cos  it  sin  (o  — N) 
sin  (T  cos  p  =  sin  u  [cos  5  sin  «'  — sin  5  cos  i  sin  (a  — N)] 

—  cos  M  sin  5  cos  (a  — N)     ...         (7) 
cos  or  =  sin  u  [sin  5  sin  t  +  cos  5  cos  i  sin  (a— N)] 

+  cos  u  cos  5  cos  (a  —  N). 

These  equations  are  perfectly  rigorous,  and  determine  o-  and  p 
without  ambiguity. 

As  the  angle  o- cannot  be  compared  directly  with  observation, 
it  will  be  more  convenient  to  express  a  in  t^rms  of  s ;  thus  from 
(6)  we  have 

(8) 


r  cos  <r  +  A 


and  since  s  is  always  very  small,  we  may  substitute  the  arc  for 
the  tangent  and  develops  the  second  member  into  a  series.  Hence 
we  have 

s- — sin  (T— i — sin2<T-i-»  — sin  3ff—     etc.  ...         (a) 

A  -  A-  ^  A'  ^^' 

Now  r  is  always  very  small  compared  with  A,  and  therefore 

the  terms  containing  the  square  and  higher  powers  of  —  may   be 

A 
neglected  without  sensibly  afiecting  the  result.     We  shall  then 
have 

/'  ^ 

*    =  : T/ "'^  ""  ••  •  ■•         (10) 


T  a" 

and  if  a"  be  the  value  of at  distance  unit  —  will  be  its 

sm  I  '  A 

value  at  distance  A ;  therefore  the  last  equation  may  be  written 

^1      1"  • 
s   =  —  sin  <T, 
A 

whence 

A 
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which,  suhstituted  in  the  first  two  equations  of  (7),  gives 

a"  (   •  \ 

s"  sin  p  = — jsin  «<  cos  t  cos  (a— N)  — cos?«  sin  (o  — N)  I         ...  ...         (12) 

s"  cos  p  =—  |sin?(  [cos  S  sin «  — sin  5 cos  i  sin  (a  — N)]  — cos  2«sin  Sees  (a  — N)  ^ 

The  error  committed  in  s"  by  omitting   the  second  term  of 

(9)  is 

,    r^    sin  2(T 

*  A'     sin  i"' 
which  is  a  maximum  when 

<r  =  45°  or  135°. 
Put 


'  A-     sin  1'' 


=  v; 


then 


T  =  n/t?  sin  2"  ; 


but  —  is  the  tangent  of  the  angle  of  greatest  elongation  as  seen 
from  the  earth ;  hence,  putting 

—  =  tan  6, 

A 

we  have 
when 


1)  = 

tan' 
sin 

2"' 

e  =  0  50, 

71  =  o-oo6 

«  =  I    50, 

7}    =    0029 

«  =  3  35. 

7J    =    01  12 

«  =  5  46, 

7J    =    0-290 

«  =  9  45. 

V  =  0-829 

Thus  we  see  that  for  the  satellites  of  Mars,  Uranus,  Nephme, 
the  three  inner  satellites  of  Jupiter,  and  the  five  inner  satellites 
of  Saturn  the  error  committed  in  s"  is  less  than  o"'3,  and  for 
Callisto,  the  fourth  satellite  of  Jupiter,  which  has  the  largest 
angle  of  elongation,  the  error  can  never  exceed  o""83.  In  order 
to  adopt  (7)  and  (12)  to  logarithmic  computation  put 

sin /cos  F  =  cos  i  cos  (a  — N) 

sin /sin  F  = —sin  (a  — N)  ...  ...         ...         (13) 

cos  /  =  —  sin  i  cos  (o  — N) 

sin  g  cos  G  =  cos  S  sin  i— sin  8  cos  i  sin  (a  — N) 

sin  y  sin  G  =  —  sin  5  cos  (o— N)  ...         ...         (14) 

COS  ff  —  cos  5  coa  i  +  siu  8  sin  i  sin  (a  —  K) 
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sin  A  cos  H  =  sin  5  sin  i  +  cos  5  cos  J  sin  (o  — N) 

sin  A  sin  H  =  cos  5  cos  (a— N)    ...         ...         ...         (15) 

cos  A  =  sin  5  cos  i  — cos  5  sin  i  sin  (a— N). 

Hence  (7)  and  (12)  become 

sin  <T  sin  p  =  sin  /  sin  (F  +  u\ 

sin  o-  cos  p  =  sia  ff  sin  (G  +  u)  16) 

cos  <r  =  sin  A  sin  (H  +  «) 

and 

s"  sin  p  =  —  sin  /  sin  (F  +  k) 
A 

s''cosp= — sin^sin(G  +  K)  ...  ...         ^17) 

A 

The  last  group,  which  is  known  as  Bessel's formula?,  gives  s' 
and  p  for  any  given  date,  and  also  enables  us  to  compare  the 
computed  place  with  that  obtained  by  observation.  When  5  and 
p  are  to  be  computed  for  a  series  of  dates  it  will  be  most  con- 
venient to  compute  the  values  of  the  auxiliaiy  angles,  /,  F, 
g,  and  G,  for  the  mean  noon  of  several  consecutive  days,  and 
then  interpolate  for  any  intermediate  date.  The  auxiliaries 
h  and  H  will  not  be  required  unless  we  also  want  the  value  of  a. 

The  following  relations  among  the  auxiliary  angles  may  be 
here  noticed,  since  they  will  be  required  presently  in  the  reduc- 
tion of  the  differential  equations. 

By  combining  the  values  of  sin  /  sin  F,  sin  g  cos  G,  sin  /' 
cos  F,  and  sin  g  sin  G,  we  have 

sin  /  sin  ^  sin  F  cos  G  =  sin*  (a  — N)  cos  i  sin  5  — sin  (a— N)  sin  t  cos  5 
sin  /  sin  g  cos  F  sin  G  =  —cos-  (a— X)  cos  1  sin  5; 

therefore  by  subtraction  we  have 

sin  /  sin^  sin  (F  — G)  =  cos  i  sin  5  — sin  (o  — N)  sin  i  cos  5 

=  cos  A       (18) 

In  a  similar  manner  we  find 

sin ^  sin  A  sin  (G  —  H)  =  cos  /      ...         ...         ...  ...         (19) 

sin /sin  A  sin  (H— F)  =  cos^r      ...         ...         ...         ...        (20) 

Again,  multiplying  the  first  and  second  of  (15)  by  cos  u  and 
sin  u  respectively,  we  have,  after  obvious  reductions — 

sin  A  cos  (H  +  u) 

=  [sin  i  sin  5  +  cos  i  co«  5  sin  (a  -  X")]  cos  u  —  cos  5  cos  (a  —  X)  sin  u  ; 

also,  multiplying  the  second  members  of  the  first  and  second 
of  (16)  by  the  third  of  (14)  and  (13)  respectively,  we  have,  by 
the  aid  of  (7) — 

cos  /  sin  g  sin  (G  +  ?<)  — sin  /  cos  g  sin  (F  +  u)  =  [sin  1  rin  8 

-r  cos  i  cos  5  sin  (o— N)]  cos  «— eos  5  cos  (o— N)  sin  u. 
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Comparing  this  with  the  preceding  equation,  we  have  the  follow- 
ing relation  : 

sin /i  cos  (H  +  m)  =  cos  /  sin^  sin  (Gr  +  «)  — sin  /  cos  ^r  sin  (F  +  m)     (21) 

To  Correct  the  Elements. 

When  the  comptited  places  do  not  agree  with  those  obtained 
by  observation  it  Avill  be  necessary  to  correct  the  elements  by 
forming  equations  of  condition,  which  will  be  of  the  form 

clu  aN  di 

dif  .         ds  .-HT     ds  ^  •      ds  . 
As  =  — Am  +  — AN  +  — Az+  — Aa, 
au  dN  di  da 

where  we  consider  for  the  present  the  elements  of  a  circular 
orbit. 

Differentiating  the  first  of  (7)  we  get  (a  and  8  being  con- 
stants) 

Bin  (T  cos^^  +  cos  a-  sinpda  =  [sin  (a  — N)  sin  u  +  cos  (o  — N)  cos  i  cos  u]du 

+  [cos  (a  — N)  cos  u  +  sin  (a  — N)  cos  i  sin  uldN 
—  cos  (a  — N)  sin  i  sin  udi 
=     sin  /  cos  (F  +  w)rfM 
+  [cos  (a—N)  cos  u 

+  sin  («— N)  cos  i  sin  u]dN     {a) 
+  cos  /  sin  udi 

Differentiating  the  second  of  (7)  gives 

—  siu  (T  sin^cZ/j  +  cos  <r  eospdir  =  [cos  (a  — N)  sin  S  sin  7«  — sin  (a  -N) 

cos  i  sin  5  cos  u  +  sin  i  cos  5  cos  7i]du 
+  [co3  (a— N)  cos  i  sin  5  sin  u 

—  sin  (a— N)  sin  S  cos  «]<ZN 
+  [sin  (a  — N)sin  i  sin  5  sin  u 

+  cos  i  cos  5  sin  n]di 
=     sin  ff  cos  {(jr  + u)d7i 

+  sin  S  sin  /  sin  (F  +  ?«)<ZN         ...         (6) 
+  cos  ^  sin  udi 

Multiplying  (a)  by  the  second  of  (16),  and  (h)  by  the  first 
of  (16).  and  subtracting  the  latter  result  from  the  former,  we 
have,  after  reducing  by  the  aid  of  the  relations  (18-21) — ■ 

sin"  ffdjo  =  —cos  //du 

+  [sin  i  sin  k  sin  (H  +  u)  —  sin  5]rZN  . . .         (22) 

+  sin  1/  cos  (II  +  11)  sin  udi. 
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Differentiating  the  third  of  (7),   we  find  after  similar  re- 
dactions 

sin  ada  =  —sin  k  cos  (TL-i- ti)du 

+  COS  5  sin  /  sin  (F+ K)rfN    ...         ...         (23) 

—cos  k  sin  udi. 

From  the  plane  triangle  EPS  (fig.  2)  we  have 

A'  =  A  cos  s  +  r  cos  (<r— s) 

and 

r     .     ,         , 
sin  s  =  —  sm  (ff—s). 

A  ' 

Differentiating  the  latter  and  redncing  by  the  former  we  have 

— — ds  =  cos  {(r  —  s)da  + sin  {a  — s)  —      ...         ...    '     (24) 

sin  s  r 

whence 

,  sin  (T  J      .       /         ^  da 

a<T  =— rfs  — tan  («r  — s)  — , 

sin  s  cos  (ff  —  s)  a 

where  we  write  a  for  r  in  the  circular  orbit. 

This   being   substituted   in    (23)    gives,    after   patting    for 
brevity 

sin  s  cos  (f  — s)  =  n, 

sin-  ads  =  —n  sin  h  cos  (H  +  ii)du 

+  B  cos  S  sin  /  sin  (F  +  u)d^ 
—n  cos  h  sin  udi 

+  - sin  <r  tan  ((r—s)<Za  (25) 

Therefore   the  partial  differential  coefficients  of  h   and  s  with 
respect  to  the  elements  of  the  circular  orbit  are 


dp 
du 

^  cr-^  =  sin  i  sin  h  sin  (H  +  «)  sin  «  — sin  5 


sin*  a-J-  =  —  cos  h 
du 


""-dN 


(26) 


sin-  <^-f.  =  sin  A  cos  (H  +  «)  sin  u. 


and 


■     ■>       ds  ■        I  -rr  V 

Sin-  a—-  =  —n  sin  h  cos  (H  +  u) 
du 

ds 

sm*  "-—  =  n  cos  5  sin  /  sin  (F  +  «) 

siu^  (T-—  =  —n  cos  n  sm  u        ...         ...         ...         ...         (27) 

-  ,     rfs       n    .  , 

8in-  ff— -  =  —  sin  (T  tan  (tr—s), 
da       a 
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where  s  and.  a  are  to  be  expressed  in  seconds  of  arc  at  dis- 
tance A. 

If,  however,  we  take  the  eccentricity  of  the  orbit  iato  con- 
sideration, "we  must  express  du  and  da  or  dr  in  terms  of  the  re- 
maining elements  of  the  elliptic  orbit.     Let 

?«o  denote  the  mean  anomaly  at  any  given  epoch  Tq, 
y.  the  mean  motion  in  a  solar  day, 
<p  the  angle  of  eccentricity  such  that  sin  <p  =  e, 
E  the  eccentric  anomaly, 
V  the  true  anomaly, 

w  the  angular  distance  from  the  ascending  node  to  the  pcriccntFe, 
and 

m  the  mean  anomaly  for  any  subsequent  date  t ; 

then  by  the  theory  of  elliptic  motion  we  have 

m  =  ?«o  +  (;!-To)/i 
'=  E  —  sin  <p  sin  E 


E 

==  m  +  (sin  (^  — 1  sin' (^  +    .  . 
+  A  sin''  d)  sin  2m—   .  .  . 

.  )  sin  m 

cos  E 

.  .  .)  cos  m 

+ 1  sin  0  cos  2m 

- 

and  also 


r  =  a(i—  sin  <p  cos  E) 

=  a  (i  —  sin  (^  cos  ?n  + ^  sin''' (^—   .  .  .)  ...         (28) 

V  =  m  +  2  sin  (p  sin  m  +  ^  sin^  (^  sin  2m  +  ... 

=  W  + w  +  2  sin  (/)  sin  m+ ...         (29) 


Since  the  eccentricity  of  the  orbits  of  all  the  satellites  of  the 
superior  planets  is  not  only  small  but  also  very  uncertain,  we 
may  neglect  the  square  and  higher  powers  of  sin  <^  without 
sensibly  affecting  the  accuracy  of  our  results. 

Difiierentiating  (29)  gives 

du  =  duo  +  dm  +  2  am  <p  cos  mdm  +  2  sin  m  cos  <\>d<p  ; 

and  expressing  the  radius  vector  of  the  satellite  in  terms  of  the 
semi-major  axis  (a)  of  its  orbit  by  putting 


we  shall  have  by  differentiating  this  last  equation 

dr       da     dp 


/ 
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Inhere  r  and  o  are  expressed  in  the  same  unit  as  A,  Therefore 
from  (28)  we  have 

dp  =  d  (- )  =  sin  (p  sin  mdm— cos  m  cos  <pd<p; 

hence 

dr  _da    sin  ^  sin  ^j„_  cos  m  cos  ^  j^ 

rap  p 

dr 
Therefore  by  substituting  these  values  of  du  and  —   in    (22), 

(24),  and  (25)  we  have  for  the  partial  differential  coeflBcients 
of  p  and  s  with  respect  to  the  elements  of  the  elliptic  orbit  the 
following  results : — 

sin*  a—  =  —cos  h 
dw 

Rin*  <r—  =  —  (i  +  2  sin  (f)  cos  m)  cos  h 
dm 

sin'  a~  =  —2  sin  m  cos  <p  coe  /*        ...         ...  ...         (30) 

dip 

sin*  (T— ^  =  sin  t  sin  h  sin  (H  +  u)  sin  «  —  sin  5 

sin*  (T— 2  =  sin  A  cos  (H  +  «)  sin  u 
di 

and 

sin*  (T— ^  =  —  n  sin  A  cos  (H  -f  k) 

sin*  a —  =  —  (i  +  2  sin  <p  cos  m)  n  sin  h  cos  (H  +  u) 
dm 

+  -  sin  <7  tan  (ff  — s)  sm  <J)  sin  w 

sin*  a —  =  —  Tzn  sin  m  sin  A  cos  (fl  +  it) 
d<p 

+  -  sin  a  tan  (<r  — 5)  cos  m]  cos  <p 

sin*  <r-~"   =  »  cos  8  sin  /  Bin  (F  +  «)  ...         ...         ...         ^31) 

.  ,     ds  ,    ■ 

sin*  «r —  =  —  »  cos  A  sin  u 
at 

•  ,     ds       n    ■        ^      ,        . 
sin*  a—-  =  -  sin  a  tan  (<r  — »), 
(2a       a 

in  which  a  is  to  be  expressed  in  seconds  of  arc  at  distance  A. 

To  Comptite  the  Times  of  greatest  Eastern  and  Western  Elongation. 

Let  us  conceive  three  planes  to  be  passed  as  follows  : — 
The  first  through  the  centre  of  the  planet  and  the  poles  of 
the  equator  and  the  orbit  of  the  satellite  ; 
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The  second  through  the  centre  of  the  planet,  the  pole  of  the 
orbit  of  the  satellite,  and  the  centre  of  the  Earth  ; 

The  third  through  the  centre  of  the  planet,  the  pole  of  the 
equator,  and  the  centre  of  the  Earth. 

These  three  planes  will  determine '  on  the  celestial  sphere  a 
spherical  triangle  having  its  centre  at  the  centre  of  the  planet. 


Let  V  A  represent  an  arc  of  the  equator,  N  S  an  arc  of  the 
orbit  of  the  satellite,  Q  the  pole  of  the  equator,  K  the  pole  of  the 
orbit  of  the  satellite,  IST  the  position  of  the  ascending  node,  P  the 
planet,  E  the  Earth. 

Join  P  E,  P  K,  P  Q,  and  in  the  plane  K  P  E  draw  P  B  per- 
pendicular to  K  P,  and  P  C  perpendicular  to  P  E  ;  then  P  B  will 
lie  in  the  plane  of  the  orbit,  and  P  C  in  the  plane  of  the  apparent 
orbit,  which  is  the  projection  of  the  real  orbit  on  the  plane  pass- 
ing through  P,  and  perpendicular  to  the  line  P  E  joining  the 
Earth  and  planet. 

Put 

e  =  the  angle  EPB, 

q  =  the  angle  of  position  of  the  minor  axis  of  the  apparent  orbit, 

Po  =  q  ±  90°,  the  angle  of  position  of  the  major  axis, 

Mq  =  the  angular  distance  of  the  end  of  the  major  axis  from  the  ascend- 
ing node  measured  on  the  orbit  of  the  satellite, 

a  and  b  =  the  semi  major  and  minor  axes  respectively  of  the  apparent  orbit ; 
then  in  the  spherical  triangle  E  K  Q  we  shall  have 

KQ  =  i,  the  inclination  of  the  orbit  of  the  satellite, 

EQ  =  90°  +  5,  the  polar  distance  of  tlie  Earth  as  seen  from  thr  planet,  and 

EQK  =  the  difference  between  the  Eight  Ascension  of  the  Earth  and  tho 
Right  Ascension  of  the  pole  of  the  orljit  of  the  satellite 

=  (a- 180°) -(N -90°) 

=  o-N-90°; 

hence  the  remaining  three  parts  of  tho  triangle  can  be  found. 
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Since  the  planes  P  K  E,  P  K  Q  pass  through  the  pole  of  the 
orbit,  they  are  perpendicular  to  its  plane,  and  the  former  inter- 
sects the  orbit  in  those  points  which,  as  seen  from  the  Earth,  are 
nearest  the  planet,  or,  in  other  words,  the  intersection  of  the 
plane  P  K  E  with  the  apparent  orbit,  determining  the  minor 
axis,  while  the  plane  P  K  Q,  being  perpendicular  to  the  planes  of 
the  equator  and  of  the  orbit  of  the  satellite,  intersects  the  latter 
at  a  point  90"  distant  from  the  node.  The  angle  E  K  Q,  being  the 
inclination  of  these  two  planes,  is  measured  by  the  arc  of  the 
orbit  intercepted  between  these  two  points  of  interaection. 

The  side  K  Q  will  intersect  the  orbit  in  a  point  whose  position 
is 

90°  +  N ; 

and  if  we  denote  the  position  of  the  intersection  of  the  side  K  E 
with  the  orbit  by  3  we  shall  have 

angle  EKQ  =  go^'  +  N-  ;8. 

But  since  the  plane  P  K  E  intersects  the  apparent  orbit  in  its 
minor  axis,  the  position  of  that  point  of  the  orbit  at  which  the 
satellite,  as  seen  from  the  Earth,  will  be  at  its  greatest  eastern 
elongation  will  be 

^3-90", 

and  therefore  the  distance  of  this  point  from  the  node  will  be 

^-9G°-N. 

Hence  we  have 

«„  =  /3-9o''-N. 

But  we  have  shown  that 

angle  EKQ  =- /3  +  90"^  +  N ; 

therefore  the  angle  E  K  Q=  —  ?/o. 

The  angle  K  E  Q  is  the  angle  of  position  of  the  minor  axis 
of  the  apparent  orbit ;  hence 

aogle  KEQ  =  q  =  p^  —  90°. 

The  side  K  E  is  the  angular  distance  of  the  Earth  from  the  pole 
of  the  orbit,  as  seen  from  the  planet,  and  is  equal  to  the  ano^le 
B  P  C  ;  therefore  we  have 

KE  =  90-  —  d, 

where   0   is   always   to  be    taken    between     +90°    and   go°. 

When  0  is  positive  the  Earth  will  be  on  the  north  side  of  the 
plane  of  the  orbit  and  the  position  angles  will  increase,  and  when 
6  is  negative  the  Earth  will  be  on  the  south  side  and  the 
position  angles  will  decrease. 
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The  fundamental  formulae  of  spherical   trigonometry  being 
applied  to  the  triangle  E  K  Q,  we  have 

cos  EK  =  cos  KQ  cos  EQ,  +  sin  KQ  sin  EQ,  cos  KQ,E 
sin  EK  cos  EKQ  =  sin  KQ  cos  EQ-cos  KQ,  sin  EQ  cos  KQE 
sin  EK  sin  EKQ  =  sin  EQ  sin  KQE 

sin  EK  cos  KEQ  =  sin  EQ  cos  KQ-cos  EQ  sin  KQ  cos  KQE 
sin  EK  sin  KEQ  =  sin  KQ  sin  KQE. 

But  we  have  shown  that 
KQ  =  i, 
EQ  =  90°  +  5, 
KQE  =  a  —  N  —  90°,  which  are  known  ; 
and 

EK  =  90° -e, 

EKQ  =  -u^, 

KEQ  =  Po  —  90°,  which  are  to  be  found. 

Substituting  these  values  of  the  sides  and  angles  in  the  pre- 
ceding formula?,  we  have 

sin  0  =  —cos  '  sin  5  +  sin  e  cos  5  sin  (a  — N) 

cos  6  cos  ?«Q  =  —sin  i  sin  5  — cos  i  cos  S  sin  (a  — N) 

cos  0  sin  Uq  =  cos  8  cos  (a  — N)      ...         (32) 

cos  0  sin  p^  =  cos  i  cos  5  +  sin  i  sin  S  sin  (a  — N) 

cos  0  cos  Po  =  sin  i  cos  (o— N) 

for  the  determination  of  0,  «„,  and  p^.     The  agreement  of  the 
two  values  of  cos  6  and  the  one  value  of  sin  d  will  serve  to  check 
the  accuracy  of  the  computation.     It  now  remains  to  transform 
these  equations  for  logarithmic  calculation. 
If  we  put 

c  sin  C  =  cos  t 

c  cos  C  =  sin  i  sin  (o  — N), 

c'  sin  C  =  sin  i 
'c'  cos  C  =  —cos  i  sin  (o  — N), 

and  remembering  that  for  the  greatest  western  elongation  —  tt^ 
becomes  —  m^+i8o°,  and  p^  becomes  ^0  + 180°,  we  shall  have 

ein  0  =  c  cos  (C  +  5) 
cos  e  cos  Mq  =  ±</  cos  (C  +  5) 

cos  0  sin  Mq  =   ±008  5  cos  (a  — N)       (33) 

cos  0  sinpo  =  ±c  sin  (C  +  5) 
cos  0  cos  Po  =   ±  Bin  i  cos  (a  -  N), 


Jaue  1884.  Satellite  of  a  Superior  Planet.  407 

in  which  the  Yipper  sign  is  to  be  used  for  the  greatest  eastern, 
and  the  lower  sign  for  the  greatest  western  elongation. 

Now  if  u  denote  the  distance  of  the  satellite  from  the  ascend- 
ing node  at  the  date  T,  and  if  t  denote  the  time  of  greatest 
elongation,  we  shall  have 


whence 

We  also  have 


(^ 

-T) 

/A 

=  Mo-W. 

t 

=  T  +  "°- 

It 

M 

b  = 

a 

COS  KE 

= 

a 

sin  e     ... 

(34) 


(35) 

which  determines  the  semi-minor  axis  of  the  apparent  orbit. 


Spottitteooile  A  Co.,  Piintrrs,  X^ir -street  S({ua/e,  London. 
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THE  ORBIT  OF  BARNAED'S  COMET,  1884. 


BY 


J.  MORRISON,   M.D.,  Ph.D. 
Assistant  on  the  American  Ephemeris,  Washington,  D.G, 


The  Orlit  of  Barnard's  Comet,  1884.  By  J.  Morrison,  M.D.,  Ph.D., 
Assistant  on  the  American  Ephemeris,  Washington,  D.C. 

From  the  nnmeroas  observations  of  this  comet,  already 
published,  I  select  for  the  purpose  of  computing  elliptic  ele- 
ments the  following  observations  made  at  Cambridge,  Mass., 
Nice,  and  Washington  respectively  : — 


Washington  If.T. 

July  25371990 
Aug.  24- 1 24442 
Sept.  23339919 


242  4  45 
267  10  423 
297  55     93 


-37  15  320 
-35  15  304 
-26  49  506 


From   these  positions   I   obtain  by   the  usual  methods   the 
following  system  of  elliptic  elements  : — 


Epoch 

M 

tt 

Q 

i 

<P 
log} 
logo 

log/1 

p 


1884.  Sept.  245  Wash.  M.T, 

7  13  »9''52\ 

300  57  4443 

^  .-,  .,^.r/;'  Mean  Equinox 

5  "  2356-     of  18840. 

5  27  1894 

35  37     250' 

o- 1069968 

0-4862043 

2  8207001 

1958-41  days. 


The  residuals  for  the  middle  place  are  zero,  and  for  the  fol- 
lowing dates  the  difference  between  the  computed  and  observed 
places  is  as  small  as  can  be  expected,  considering  the  great 
difficulty  in  observing  accurately  so  faint  an  object : — 


Wash. 

M.T. 

dXooep. 

c-o. 
dp. 

log.  A. 

1884,  Aug. 

12-361652 

H 

+    6  52 

-  i'zi 

965375 

Sept. 

15-408454 

+    985 

-  809 

974018 

Oct. 

II-3I8804 

+  30-9* 

— 10-62 

986241 

14-365442 

+  21-56 

+  15-71 

987849 

^ottweoode  i  Co.  Printert,  Iffvc-ttntt  Square,  London. 
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THE  ORBITS   OF   COMETS  FABRY  AND  BARNARD- 
HARTWIG. 


EPHEMERIDES  OF  THE    SATELLITES  OF  MARS 
DURING  THE  OPPOSITIONS  OF  1888  AND  1890. 

BY 

J.  MORRISON,  M.D.,  Ph.D. 

Assistant  on  the  American  Ephemeris,  and  Professor  of  Chemistry 
in  the  National  University,  Washington. 
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The  Orbits  of  Comets  Fabry  arid  Barnard- Hariwig.  By  J. 
Morrison,  M.D.,  Ph.D.,  Assistant  on  the  American  Ephemeris, 
and  Professor  of  Chemistry  in  the  National  University, 
"Washington. 

Comet  Fabry. 

The  observations  npon  vrhich  the  following  hjperbolic  ele- 
ments of  this  comet  are  founded  are  as  follow  : — 


Greenwich  M.T. 

Appar.  a 

Appar.  S 

d 
1885,  Dec     7536032 

h   m       s 
0  24  460 1 

+  20  52' 34-9 

1886,  Mar.    7"3i6539 

23   19  4292 

+  31    16  44-8* 

„     April   1-339713 

23   18  3086 

+  38  37  570 

„     June    6-953118 

8  47  5215 

-40  38    0-5 

The  first  is  a  meridian  observation  made  at  Ann  Arbor ; 
the  second  results  from  extra  meridian  observations  made  at 
Greenwich  and  Paris,  the  comparison  star  being  the  same  at 
both  places ;  the  third  was  made  at  Bothkamp  and  was  obtained 
from  Astronomiscite  Nachrichten,  Xo.  2703;  and  the  fourth  was 
made  at  Sydney,  the  comet  being  at  the  time  "  extremely  faint, 
but  in  a  good  position  for  observation  with  Cape  Cat.  (1880) 
4707,"  fifty  comparisons  having  been  made.  (Monthly  Noticet; 
vol.  xlvi.  p.  496.)  These  observations  were  corrected  for  aber- 
ration and  pai"allax  by  means  of  approximate  parabolic  elements. 

*  Geocentric. 
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T 

d 
1886,  April  5-9520  Greenwich  M.T. 

«ii 

0      (       // 
126  34  49215 

Q, 

36  22  11-454  Mean  Equinox,  1886 

i 

82  37     6-012 

e 

1-00047857 

log  a 

3-1278354 

log  q 

98077809 

The  two  middle  places  are  well  represented.* 
The   interval    between   the  extreme   observations  is   181-42 
days,  during  which  the  comet  described  208°  9'  of  its  orbit. 


Comet  Barnard- Hartwig.     1886. 

From  the  following  observations  the  first  and  third  of  which 
were  made  at  Washington,  D.C.,  and  the  second  at  Kiel  {Astj-ono- 
mische  Nachrichten,  2753),  an  elliptic  orbit  of  long  period  is 
obtained,  the  elements  of  which  are  given  below  : — 

Greenwich  M.T.  Appar.  o.  Appar.  S. 

d  h     m        s  o       ,       II 

1886,  Oct.    7-918987  10  42    9-32  +    I  22    63 

„       „    29708660  II  39  2282  +   5  49  I9"4 

„    Dec.   2-976975  15  32    391  + 17  58  55-8 

By  means  of  approximate  parabolic  elements  computed  from 
a  shorter  interval,  the  corrections  for  aberration  and  parallax 
were  obtained  and  applied. 

d 
T  1886,  Dec.  16-514158  Greenwich  M.T, 

O  I  II 

w  86  21  58- 570 

Q,         137  21   36- 163  Mean  Equinox,  18860 
i  78   22   25-525 

e  0-9987252T 

log  a  2-7162x51 

log }  98216538 

Motion  retrograde.  The  middle  place  is  exactly  represented. 
These  elements  give  a  period  of  11,866  years,  which  is  of  course 
very  uncertain,  since  the  interval  between  the  extreme  observa- 
tions is  far  too  short  to  determine  this  element  accurately  in  an 
orbit  such  as  this  is. 

The  formulae  for  the  equatorial  rectangular  coordinates  are — 

x  =  [9-8740o86]r8in(D+  6  52  45-171) 
2/  =  [9-826898o]rsin(y+i98  35  9651) 
c  =  [9  9977333]  r  sin  (».+ 102     5  28657) 
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Ephemerides  of  the  Satellites  of  Mars  during  the  Opposition.'^  of 
1888  and  1890.  By  J.  Morrison,  M.D.,  M.A.,  Ph.D., 
Assistant  on  the  American  Ephemeris,  and  Professor  of 
Chemistry,  National  University,  Washington, 

These  ephemerides  have  been  computed  from  the  following 
elements  of  the  orbits  of  the  satellites,  referred  to  the  equator 
and  equinox  of  the  respective  epochs : — - 


Phobos. 

Epochs 

1888,  April  ii-o 

1890,  May  270  Greenwich  il.T. 

Period 

03189113  (mean  solar) 

M 

1128-8405 

a 

12-953  (at  distance  unity) 

i 

36  44-6 

0       / 
36  44  0 

N 

47  i8-2 

47  19*0 

» 

17  40 

Deimos. 

117  53-6 

Epochs 

1888,  April  no                    ] 

[890,  May  270  Greenwich  M.T. 

Period 

d 

1-262435  (mean  solar) 

M 

2S5- 16322 

a 

32-354  (at  distance  unity) 

i 

35  36 

35  35-'s 

N 

48  io'66 

48  11-5 

u 

225  41 

172  21-5 

Crreenwich  Mean  Time  of  Greatest  Elongation. 
Phobos. 


G.M.T. 

P 

a 

6 

G.M.T. 

P 

a 

1888.    d 
Mar.  21 

h    m 

7  535  W 

306-4 

I9'i 

6' I 

188S.    d 
Mar.  30 

h      m 
6    85  W 

0 

II 

22 

ID  40-4  E 

31 

8  55-3  E 

23 

13  273  W 

Apr.     I 

II  421  W 

305-8 

20-5 

24 

16  14-2  E 

2 

14  290  E 

25 

19     i-i  W 

3 

17  15-8  W 

26 

21  480  E 

126-2 

19-6 

6-4 

4 

20    2-6  E 

28 

0  34 -8  W 

5 

22  49-4  W 

29 

3  21-7  E 

7 

I  36-1  E 

1253 

210 

6-8 


73 
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G.M.T.  p             a         b                               G.M.T.                 p             a          b 

1888.    d      h       m  o             II         II           1888.    d    h        m                  q             hi, 

Apr.    8    4  228  W  Apr.  24  22     4-4  E 

9  7     96  E  26    o  5i'3  "W 

10  9  563  W  27     3  382  E 

11  12  431  E  28    6  251  W 

12  IS  298  W  3049    213     76               29    9  11-9  E     1233     210     80 

13  18  166  E  30  II  588  W 

14  21     3-3  W  May    I   14  45-8  E 

15  23  50-1  E  2  17  327  W 

17  2  36-9  W  3  20  197  E 

18  5  237  E  124-4     21-4      78                 4  23    67  W    3027     205      8-1 

19  8  105  W  6     I  53-8  E 

20  10  57-2  E  74  40-9  W 

21  13  440  W  -                               87  28-0  E 

22  16  308  E  9  10  17-1  W 

23  19  17-6  W  3037     21-2     8'0               10  13    4*3  E     123-0    20-0     8-2 


Deimos. 

G.M.T.  p            a           b                           G.M.T.               p            a           b 

1888.    d      h      m  o              II          II            i888.    d     h       m                   ^             n         ,i 

Mar.  20  15  43-4  E  1251     478    15-5      Apr.  23  16  47-0  E 

22  13     7-9  W  25  14  io'7  W 

24  10  32-3  E  27  II  34-4  E     122-1     52-6    20-1 

26     7  56-5  W  •                       29     8  583  W 

28     5  20-5  E  May    i     6  22-4  E 

30    2  444  W  304-6     50-4    168                 33  469  W 

Apr.    I     o    8-1  E  51  11-9  E 

2  21  31-6  "W  6  22  37-2  "W   301-4    51-2   20-1 

4  18  54-9  E  8  20    2-8  E 

6  16  18-2  W  10  17  28-6  W 

8  13  41-4  E  124-0    52-6    18-4               12  14  54-6  E 

10  II     4-5  W  14  12  20-8  W 

12    8  27-6  E  16    9  47-1  E     121-0    50-6   200 

14    5  50-7  W  18     7  13-4  W 

16    3  13-8  E  20    4  39-8  E 

18  o  37-0  W  302-9     53-5    19-6              22    2    6-4  W 

19  22    0-3  E  23  23  33*1  E 

21  19  23-6  W  2$  21     0-0  W   300-2    48 8    197 
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Phobos. 

G.iLT. 

J> 

a 

b 

GJLT. 

p 

a 

h 

1890.     d 
May    8 

h     m 
0  448  E 

I23°6 

23-5 

2-3 

1890.    d 
June    7 

h      m 

3  521  W 

- 

" 

M 

9 

3  318  W 

8 

6  39-0  E 

10 

6  188  E 

9 

9  260 W 

II 

9     5-8  W 

10 

12  12-9  E 

1267 

265 

53 

12 

II  528  E 

II 

14  599  W 

13 

14  39  8  W 

304-1 

25-0 

30 

12 

17  468  E 

14 

17  268  E 

13 

20  338  W 

15 

20  138  W 

14  23  208  E 

16 

23    08  E 

16 

2    77  W 

3071 

260 

5-6 

18 

1  478  w 

17 

4  54-8  E 

19 

4  347  E 

1245 

254 

33 

18 

7  41  7  W 

20 

7  217  W 

19 

10  28-9  E 

21 

10    8  6  E 

20 

13  160  W 

22 

12  556  W 

21 

16    32  E 

1273 

25-5 

57 

23 

15  425  E 

22 

iS  50  2  W 

24 

18  295  W 

3050 

260 

3-8 

23 

21  37-4  E 

25 

21  16-4  E 

25 

0  246  W 

27 

0    3-3  W 

26 

3  11-9  E 

28 

2  502  E 

27 

5  59-2  W 

3075 

24-9 

5-7 

29 

5  370  W 

28 

8  46  4  E 

30 

8  239  E 

1256 

263 

4'3 

29 

II  33  7  W 

31 

II  IO-8  W 

30 

14  2roE 

June    I 

13  57-6  E 

July 

I 

17    8-3  W 

2 

16  446  W 

2 

19  556  E 

1277 

241 

57 

3 

19  314  E 

3 

22  42-9  W 

4 

22  184  W 

3061 

267 

4-8 

5 

I  302  E 

6 

I     5-2  E 

6 

4  17-6  W 

Deimos. 

1890.     d 
Apr.  30 

G.M.T. 
h     m 
22   40-9  E 

p 
1225 

a 
55'i 

b 

1890.    d 
May  16 

OJC.T. 
h      m 
2     36  E 

p 

a 

6 

May    2 

20     69  W 

17 

23  280  W 

4 

17  327  E 

19 

20  523  E 

1240 

63-7 

94 

6 

14  583  w 

' 

21 

18  166  W 

8 

12  236  E 

23 

15  409  E 

10 

9  489  W 

3031 

597 

7-2 

25 

13   51  w 

12 

7  140  E 

27 

10  29-1  E 

14 

4  38  9  W 

29 

7  53  I  W 

304-9 

661 

irS 
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G.M.T. 

p 

a 

6 

1890.    d 
June  19 

G.M.T. 

V 

1890.    d 
May  31 

h       m 

5  17-iE 

0 

" 

" 

h    m 
3  21-4  E 

0 

June   2 

2  411  W 

21 

0  467  W 

4 

0    53  E 

22 

22  123  E 

5 

21  29-4  W 

24 

19  381  W 

7 

18  536  E 

1256 

666 

13-8 

26 

17     4-1  E 

1265 

9 

16  179  W 

28 

14  30-3  W 

II 

13  424  E 

30 

1 1  566  E 

13 

II     69  W 

July    2 

9  23-1  W 

15 

8  31-5  E 

4 

6  498  E 

17 

5  56-3  W 

306- 1 

656 

14-8 

6 

4  i6-8  W 

3066 

621    15-2 


57-6    145 


In  these  ephemerides  ip  denotes  the  position-angle  of  the 
major  axis  of  the  satellite's  apparent  orbit,  and  is  reckoned  from 
the  north  towards  the  east,  and  from  0°  to  360°  ;  a  and  h  denote  the 
major  and  minor  semi-axes  of  the  apparent  orbit.  During  these 
two  oppositions  of  the  planet  the  satellites  move  in  the  direction 
of  increasing  position-angles,  the  Earth  being  above  the  plane  of 
the  orbits.  The  time  of  greatest  elongation  has  been  given 
to  the  nearest  tenth  of  a  minute  in  order  that  a  sufficiently  accu- 
rate comparison  may  be  had  with  observation. 


.'Spottl.vroolie  <fc  Co.  Printers,  Nfie-itreet  Square,  London. 
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